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Abstract 

It  will  be  shown  that  some  important  characteristics  of  stationary  operating  fuel  cells  can  be  unCnstood  by  analyzing  a  simple  basic 
isothermal  cell  model.  The  model  is  independent  of  the  type  of  fuel  cell  and  will  be  used  for  a  genera!  description  of  cells  supplied  with 
gaseous  fuel  and  oxidant.  An  example  is  also  given  in  which  the  results  are  compared  with  published  data  of  a  bench-scale  molten  carbonate 
fuel  cell  (MCFC).  Simple  expressions  relating  the  cell  voltage,  cell  current  and  fuel  utilization  are  derived  and  found  to  be  accurate  despite 
their  simplicity.  Deviations  in  analytically  calculated  cell  voltages  from  numerical  computations  are  in  the  order  of  mV.  Moreover,  it  will  be 
shown  that  a  non-homogeneous  current  distribution  yields  an  extra  overpotential  in  the  older  of  10  mV  for  an  MCFC. 
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1.  Introduction 

Engineering  aspects  often  dominate  the  setup  of  fuel  cell 
models  and  make  them  rather  complicated.  Therefore,  most 
fuel  cell  models  are  not  suitable  for  an  analytical  approach 
and  require  extensive  computer  software  for  generating 
numerical  results.  In  addition,  the  dependence  of  fuel  cell 
behaviour  on  the  model  parameters  is  often  difficult  to 
express  on  the  basis  of  numerical  results  only.  The  great 
advantage  of  an  analytical  model  is  that  relations  between  the 
variables  can  be  made  clear  in  a  direct  way  and  that  fuel  cell 
performance  can  be  determined  by  the  derivation  of  analytical 
expressions.  On  the  other  hand,  simple  mathematical  analyses 
arc  only  possible  for  highly  idealized  models.  Nevertheless, 
analytical  results  arc  important  in  the  explanation  of  results 
obtained  from  complicated  numerical  models  as  well  as  of 
empirical  relations  obtained  from  measurements  So,  despite 
their  limitations,  analytical  results  contribute  to  a  better 
understanding  of  fuel  ceil  behaviour. 

In  this  paper,  a  simple  basic  fuel  cell  model  is  introduced 
for  a  rectangular  flat  plate  fuel  cell.  The  objective  of  the  'basic 
model’  is  to  provide  insight  in  the  elementary  characteristics 
of  stationary  operating,  isothermal  fuel  cells,  using  simple 
mathematics.  In  the  ‘basic  model’  a  fuel  cell  is  described  as 
an  equivalent  electrical  circuit  (see  Fig.  1 ).  This  approach  is 
also  used  in  Refs.  [  1—3]  and  will  be  elaborated  into  more 
detail  in  this  paper.  Special  attention  will  be  paid  to  a  linear 
relation  between  the  cell  voltage,  the  total  fuel  utilization  and 
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Fig.  1 .  Representation  of  a  fuel  cell  as  an  equivalent  electrical  circuit. 

the  cell  current.  Although  measurements  on  a  bench-scale 
molten  carbonate  fuel  cell  (MCFC)  confirm  this  double  lin¬ 
ear  relation  [  1  ],  a  detailed  derivation  of  the  relation  is  not 
given  in  any  of  the  references.  In  this  paper,  a  detailed  deri¬ 
vation  will  be  given  and  it  will  be  shown  that  more  accurate 
results  arc  obtained  if  the  relation  is  slightly  adapted. 

Table  1  shows  the  list  of  assumptions  that  underlie  the 
setup  of  the  ‘basic  model’.  Oxidant  utilization  can  be  taken 
into  account  in  a  similar  way  as  fuel  utilization  ( see  Appendix 
A),  but  will  be  neglected  for  reasons  of  simplicity.  Due  to 
the  assumptions  5  and  6,  the  setup  of  the  ‘basic  model’  is 
possible  without  paying  special  attention  to  the  theoty  of 
chemical  kinetics.  Despite  the  assumptions  listed  in  Table  1, 
the  ‘basic  model’  can,  in  most  cases,  only  be  ‘solved'  by 
using  a  numerical  method.  However,  it  will  be  shown  that 
linearization  of  the  Nernst  potential  in  the  ‘basic  model’ 
yields  an  ‘analytical  cell  model’,  suitable  for  the  derivation 
of  analytical  expressions  that  are  in  good  agreement  with 
numerical  computations  and  measurements. 
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Table  I 

Overview  of  ihe  model  assumpuons 


t  Stationary  fuel  celt  <  time  independent ) 

2  Isothermal  fuel  cell 

3  Large  oxidant  flow  {negligible  oxidanl  utilization) 

4  Changes  in  fuel  composition  are  significant  in  the  direction  of  the 

cell  outlet  only 

5(a)  Local  overpotcntial  depends  linearly  on  the  local  current  density 
Ihrough  a  quasi-ohmic  resistance,  r 
5(b)  Quasi-ohmic  resistance,  r.  is  independent  of  position 

6  All  reactions  in  the  gas  phase  are  in  equilibrium 

7  Nemst  potential  is  considered  to  be  independent  of  hydrostatic 

pressure  gradients 


2.  Basic  isothermal  fuel  cell  model 

The  total  anode  surface  of  the  cell  is  given  by  BL,  where 
L  and  B  are,  respectively,  the  length  and  the  width  of  the  cell. 
The  approach  will  be  one-dimensional,  which  means  that  ail 
variables  are  considered  as  a  function  of  the  X-eoordinate 
only,  i.e.  the  horizontal  distance  (Fig.  1 )  between  a  cross 
section  of  the  cell  and  the  cell  inlet. 

Changes  in  fuel  and  oxidant  composition  occur  due  to  the 
electrochemical  reactions  that  are  taking  place  inside  the  ceil. 
Hence,  the  local  equilibrium  potential,  as  given  by  Ncrnst 
law,  depends  on  the  local  composition  of  fuel  and  oxidant. 
The  local  overpotential  tj(X)  is  defined  as  the  difference 
between  the  local  equilibrium  potential  (V^)  and  the  cell 
potential  ( VL.C„)  at  which  the  current  is  delivered.  Assumption 
5(a)  in  Table  1  leads  to 

r)(X)=Vtiti(X)-V„H  =  ri(X)  (I) 

This  linear  current-voltage  relation  is  not  evident  frnm  a 
theoretical  point  of  view.  The  quasi-ohmic  resistance,  r,  does 
not  only  account  for  ohmic  losses,  but  also  for  voltage  losses 
due  to  concentration  and  activation  polarization.  Eq.  ( 1 )  is 
realistic  for  an  MCFC  at  normal  operation,  because  experi¬ 
mental  current-voltage  relations  arc  almost  perfectly  linear 
for  this  type  of  cell  (sec  e.g.  1 1 1 ).  But  for  low-temperature 
fuel  cells  a  strong  decrease  in  cell  voltage  is  observed  at  low 
current  densities,  usually  ascribed  to  the  activation  polariza¬ 
tion.  Furthermore,  at  high  current  densities  diffusion  limita¬ 
tions  occur  14.5).  Nevertheless,  the  current-voltage  relation 
of  a  low-temperature  cell  can  be  linearized  accurately  for 
intermediate  current  densities,  yielding  an  expression 
for  jj(X)  as  V'^+rifX).  Rewriting  this  expression  gives 
Eq.  ( 1 )  with  an  adapted  value  for  V^X).  Hence,  for  low- 
temperature  cells,  the  initial  value  1^(0)  must  be  set  equal 
;o  the  opcn-circuil  voltage  (OCV)  minus  the  correction  term 
Vcor. 

Note  that  the  resistance  of  the  current  collectors  is  assumed 
to  be  negligible  in  the  definition  of  the  overpotential  17  in 
Eq.  ( 1 ),  so  that  the  voltage  between  two  current  collectors  is 
equal  to  V...„  throughout  the  cell.  For  reasons  of  simplicity, 
we  assume  also  that  the  quasi-ohmic  resistance  is  uniform, 
i.c.  r(X)  -r  for  all  X.  In  most  cases,  the  quasi  ohmic  resis¬ 
tance  depends  on  temperature  mainly  and,  hence,  this  sim¬ 


plification  is  a  direct  consequence  of  the  assumption  that  the 
fuel  cell  is  isothermal.  In  general,  electrode  polarization 
depends  also  on  the  local  gas  compositions,  thus  on  X.  But 
in  this  paper  we  will  also  neglect  this  influence,  which  is  often 
of  secondary  importance. 

Below  we  will  introduce  an  ordinary  differential  equation 
( ODE)  for  the  fuel  utilization,  u,  based  on  the  expression  for 
the  overpotential  given  in  Eq.  ( 1 ).  The  equilibrium  potential 
will  appear  in  this  ODE  as  a  function  of  u  and  not  as  a 
function  of  X.  In  fact,  the  equilibrium  potential  depends 
mainly  on  the  local  fuel  and  oxidanl  composition  and  is  prac¬ 
tically  independent  of  the  small  total  pressure  gradients  in  the 
process  flows  [4].  Neglecting  hydrostatic  pressure  gradients 
it  is  possible  to  determine  the  composition  of  the  process 
flows  as  a  function  of  the  fuel  utilization,  it.  The  Nemst 
equation  gives  thus  the  local  equilibrium  potential  in  terms 
of  the  utilization  u(X). 

For  reasons  of  simplicity,  wc  assume  a  large  flow  of  oxi¬ 
dant  being  considered  as  homogeneous.  Thus,  we  have  to 
take  into  account  the  utilization  of  fuel  only.  It  is  convenient 
to  represent  the  maximum  amount  of  current,  that  is  produced 
in  the  hypothetical  situation  of  complete  fuel  utilization,  by 
which  we  will  call  the  ‘equivalent  input  current’.  Note 
that  /,„  can  be  calculated  from  the  gas  composition  and  the 
gas  flow  rate,  using  Faraday’s  law.  The  delivered  current  is 
denoted  by  7oul  and  will  be  smaller  than  the  equivalent  input 
current  /,„.  The  ratio  /oul//,n  gives  the  total  fuel  utilization 
and,  in  an  analogous  way,  the  local  fuel  utilization  «(X)  is 
defined  as  the  fraction  of  /,„  generated  in  the  area  between  the 
cell  inlet  and  position  X 
x 

u(X)  d£,  (0£m£1)  (2) 

where  B  is  the  width  of  the  cell.  The  expression  for  the 
overpotential  in  Eq.  ( 1 )  can  be  used  to  eliminate  the  local 
current  density  r  in  Eq.  (2).  Differentiating  the  result  with 
respect  to  X  and  writing  V^OO  instead  of  V^fX)  we  obtain 
the  following  ordinary  differential  equation 

du  B 

^=-T(V04(h)-K-ci1),  (0  <X<L)  (3) 

dX  rlm 

From  the  definition  of  the  utilization  u  in  Eq.  (2)  it  follows 
that  u  is  zero  at  the  cell  inlet,  i.e.  for  X=0.  Apart  from  this 
initial  value  for  w,  the  ODE  must  be  supplied  with  either  a 
suitable  end  value  for  u  or  a  value  for  V,*i-  Choosing  the  first 
possibility,  i  t  imposing  an  end  value  u,  for  the  utilization, 
we  have  to  solve  Eq.  (3)  for  u(X)  and  VccU,  taking  into 
account  the  boundary  conditions 

«(0)  =0,  «(/,)=«„  (0<u,<l)  (4) 

In  order  to  facilitate  the  evaluation  of  the  results,  it  is  con¬ 
venient  to  scale  Eqs.  (3)  and  (4).  To  this  end,  we  succes¬ 
sively  relate  the  equivalent  input  current  /,„  to  the  cell  surface 
BL  and  the  space  coordinate  X  to  the  cell  length  L.  The  scaled 
variables  are  denoted  by  i,n  and  x  respectively 
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=  Jr=7.  (0<jc<  l)  (5) 

DL  Li 

Scaling  Eqs.  (3)  and  (4),  by  using  these  scaled  variables, 
we  obtain  the  following  boundary  value  problem  ( B VP) 

Ih(0)=0,  U(l)  =K,  ^Anu/fin 

In  this  BVP  the  first-order  differential  equation  for  the  utili¬ 
zation  is  supplied  with  two  boundary  conditions.  This  implies 
that  the  solution  of  the  BVP  is  twofold,  i.e.  is  given  by  the 
set  («(*),  Vccl, ) .  This  solution  depends  n  the  product  riin 
rather  than  on  the  individual  data  for  r  and  iin.  However,  the 
local  current  density  does  depend  on  these  individual  data, 
e.g.  the  average  current  density  over  the  whole  cell  is  equal 
to 

Once  the  BVP  has  been  solved,  the  local  current  density  i 
can  be  calculated  from 


u.  Oxidant  utilization  in  an  MCFC  is  often  sufficient  low  to 
have  little  overall  effect  [  6] ,  because  in  practice  large  oxidant 
flows  are  used  for  cooling.  Therefore  and  for  reasons  of  sim¬ 
plicity,  it  is  convenient  to  neglect  oxidant  utilization  in  the 
‘basic  MCFC  model’.  The  representation  of  the  equilibrium 
potential  inside  an  MCFC  for  zero  oxidant  utilization 

is  derived  in  Appendix  A,  assuming  temperature  and  pressure 
to  be  uniform.  Furthermore,  we  do  not  consider  the  methane 
reforming  reaction  and  assume  that  the  water-gas  shift  reac¬ 
tion  is  in  equilibrium.  The  expression  resulting  for  V^fii) 
turns  out  to  be  complicated  and  unsuitable  for  analytical  pur¬ 
poses;  it  can  be  used  in  numerical  computations,  whose 
results  are  given  in  the  next  sections.  In  order  to  make  a  simple 
mathematical  analysis  of  the  ‘basic  model*  possible,  a  linear 
approximation  for  the  equilibrium  potential  will  be  derived 
below. 

The  Nernst  law  gives  the  equilibrium  potential  in  terms  of 
partial  pressures  p,  of  the  species  in  the  process  flows 


=£U+§,n(**3'Co|) +f 1 


(7) 


(0<x<  1)  (6) 

dx 

Note  that  the  individual  values  for  1^(0)  and  Vcen  are  irrel¬ 
evant,  because  only  the  overpotential  <  V^w)  -  VtfcU) 
appears  in  the  BVP.  Hence,  the  difference  0)  —  Vccli  is 
important;  it  is  this  difference  rather  than  the  individual  terms 
that  influence  the  local  current  density  i  and  the  local  fuel 
utilization  u. 

In  this  paper,  the  ‘basic  model'  is  only  used  to  calculate 
the  electrical  variables,  i.e.  the  current  density,  i,  the  fuel 
utilization,  u( x),  and  the  cell  voltage  V^,,.  We  distinguish 
this  set  of  variables  from  the  set  of  so-called  hydrodynamical 
variables.  In  fluid  dynamics  the  term  hydrodynamical  varia¬ 
bles  refers  to  the  hydrostatic  pressure,  the  fluid  density  and 
the  velocity  of  a  fluid  flow.  In  afuel  cell  model,  the  calculation 
of  hydrodynamical  variables  requires  a  detailed  description 
of  the  configuration  of  the  gas  channels.  However,  as  long  as 
sufficient  gas  flow  is  allowed  for,  these  details  appear  to  be 
irrelevant  for  the  calculation  of  the  variables  i,  u  and  V^,,, 
within  the  ‘basic  model’. 

Solving  the  BVP  requires  an  expression  for  the  local  equi¬ 
librium  potential  in  terms  of  the  utilization  u.  To  obtain  an 
expression  for  V^a)  we  simply  assume  fuet  and  oxidant  to 
be  a  mixture  of  ideal  gases  and  neglect  concentration  gradi¬ 
ents  over  a  cross  section  of  a  gas  channel  (ideal  mixing  in 
the  gas  phase). 

3,  The  linearized  Nernst  potential  for  a  ‘basic  MCFC 
model’ 

In  this  section  the  general  ‘basic  model’  is  specified  for  an 
MCFC.  This  requires  the  representation  of  the  equilibrium 
potential  inside  an  MCFC  as  a  function  of  the  fuel  utilization 


where  £°  is  determined  by  thermodynamics  and  the  super¬ 
scripts  ‘a’  or  V  refer,  respectively,  to  the  anode  or  cathode 
gas  composition.  As  we  assume  a  homogeneous  cathode  flow, 
the  second  term  on  the  right-hand  side  of  Eq.  (7)  is  consid¬ 
ered  to  be  independent  of  position.  In  Appendix  A  it  is  shown 
■:ow  the  partial  pressures  in  the  last  term  of  Eq.  (7)  can  be 
expressed  in  terms  of  the  fuel  utilization  u. 

At  650  °C  the  factor  RT/2F  in  the  Nernst  equation  is  equal 
to  about  40  mV,  which  is  small  in  comparison  with  common 
values  for  the  OCV  (typical  1000  mV).  The  Nernst  potential 
Vej(u)  in  Eq.  (7)  varies  thus  moderately  with  «,  except  if 
any  of  the  partial  pressures  in  the  last  logarithmic  term  is  very 
low.  The  shift  reaction  will  not  cause  extreme  low  partial 
pressures  due  to  its  equilibrium  state.  So,  only  for  values  of 
h  near  zero  or  unity  we  may  expect  big  changes  in  the  loga¬ 
rithmic  term.  Nevertheless,  Fig.  2  shows  that  for  most  values 
of  u,  an  accurate  approximation  V£,(ti)  for  the  local  equilib¬ 
rium  potential  V,ei(«)  can  be  obtained  from  a  linearization 

V*(«)  =  V*(0)-a«, 
dV 

a = constant  = - -21  for  most  (0,1)  (8) 

d  u 

The  constant  of  proportion  a  in  the  linearized  equilibrium 
potential  is  not  exactly  defined,  but  must  be  a  characteristic 
value  for  the  derivative  of  —  V<M(«)  with  respect  to  u.  Note 
that  the  dimension  of  a  is  V,  since  u  is  dimensionless. 

As  shown  in  Fig.  2,  the  local  equilibrium  potential  V^k) 
inside  an  MCFC  contains  an  initial  dip  near  the  origin,  i.e. 
the  equilibrium  potential  decreases  more  rapidly  near  «  =  0 
than  in  the  neighborhood  of  «=  1/2.  Moreover,  the  equilib¬ 
rium  potential  also  decreases  fast  near  the  point  of  singularity 
h=  1.  Apparently,  a  suitable  value  for  a  is  obtained  if  the 
derivative  of  the  equilibrium  potential  is  determined  in  the 
point  of  inflection,  i.e.  the  point,  say  k*,  where  the  second 
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Fig.  2.  Example  of  a  linearized  Nernstequalion  V*(u)  in  case  of  an  MCFC 
Note  the  initial  dip  in  the  exact  Nemsl  equation  VCI|(«)  at  low  utilization. 
( - )  (V'5|(u),and  (— - )  Vn,(M):  for  parameter  values  see  Table  2. 


derivative  of  the  equilibrium  potential  with  respect  to  u  is 
zero.  Then  we  have 


d2V  dV 

0, 

dtc  dw 


V^(  0)  =  VUI|(  «* )  +  an*  (8a) 

For  instance,  consider  the  inlet  fuel  composition  presented  in 
Table  2  and  a  homogeneous  cathode  flow  of  70%  air  and  30% 
CO 2-  In  Table  2  the  mole  fractions  after  equilibrium  are  based 
on  a  value  of  K =0.5  for  the  equilibrium  constant  of  the 
reverse  shift  reaction  ( corresponding  to  a  cell  temperature  of 
650  °C  and  atmospheric  pressure  [7] ).  For  the  fuel  compo¬ 
sition  presented  in  Table  2.  the  point  of  inflection  of  the 
function  Vctl(«)  can  be  determined  numerically  to  be 
u*  =  0.50.  In  this  point  the  slope  of  VctJ(n)  is  0.18  V.  Hence 
in  this  case  the  linear  approximation  for  the  equilibrium 
potential  is  given  by 

V*(k)  =  1.03 -0.18k,  (0.2  <u  <0.8)  (9) 

where  the  actual  OCV  for  this  gas  composition  is  1.06  V. 
Because  of  the  initial  dip  in  the  Nernst  potential,  the  value  of 


V^(0)  is  30  mV  lower  than  the  actual  CCV.  The  boundaries 
for  u  given  in  Eq.  (9)  arc  not  strict,  but  give  an  idea  of  the 
interval  for  which  the  linear  approximation  is  accurate. 

If  the  local  equilibrium  potential  of  a  low-temperature  cell 
is  linearized,  then  the  difference  V£,(0)  —1^(0)  must  not 
be  confused  with  the  correction  term  V,.ir  that  was  introduced 
in  the  previous  section. 

Machiclse  [  1  ]  determined  experimentally  the  quasi -ohmic 
resistance  r  for  an  MCFC  to  be  1 . 1 95  £  cm2,  see  Table  2.  He 
also  applied  the  PSI  (Physical  Sciences  Inc.)  cell  model  [8] 
for  his  experimental  conditions  and  obtained  a  10%  lower 
value  of  1.08  ft  cm2.  Throughout  th  s  paper  we  will  use 
r—  1.08  ft  cm2.  Only  in  Section  5,  where  a  comparison  is 
made  with  the  experimental  results  obtained  by  Machiclse. 
wc  will  use  r  =  1 . 1 95  ft  cm2. 


4.  Analytical  solutions  of  the  ‘ basic  model’ 


The  BVP  that  was  introduced  in  Section  2  can  only  be 
solved  in  closed  form  for  some  simplified  versions  of  the 
equilibrium  potential  Vcij(k).  In  this  section  an  analytical 
solution  of  the  BVP  is  derived,  assuming  the  function  Vcq(  k) 
to  be  linear.  In  Section  3  it  was  shown  for  an  MCFC  that 
VC4(k)  is  almost  linear  up  to  high  values  for  the  fuel  utiliza¬ 
tion  u.  The  analytical  solution  will  be  thus  fairly  accurate  up 
to  high  values  for  the  total  fuel  utilization  m,.  The  obtained 
analytical  expressions  will  therefore  apply  to  problems  with 
practical  values  for  k,. 

Also  approximate  solutions  for  the  BVP  will  be  derived 
and  in  Section  4.2  the  influence  of  the  initial  dip  in  the  equi¬ 
librium  potential  VcqOO  will  be  considered. 

4. 1.  Analytical  cell  model 


The  cell  model  that  is  obtained  if  a  linear  function  is  used 
to  describe  the  equilibrium  potential  V^u)  in  the  BVP 
is  referred  to  as  the  ‘analytical  cell  model’.  The  problem 
is  to  find  a  function  u:  [0,1  ]  ->  [O.k,]  and  a  constant  V£„ 
satisfying 


Table  2 

Standard  gas  compositions  for  an  MCFC  and  corresponding  parameter  values  at  atmospheric  pressure;  PSI  --  calculated  from  PSI  model  [  1  ] ,  exp.  =  determined 
from  experimental  results  11].,  Methane  reforming  equilibrium  is  neglected.) 


Anode  inlet  gas 

H, 

CO 

H,0 

CO, 

ch4 

Before  equilibrium 

0.64 

0 

0.20 

0.16 

0 

After  equilibrium 

0.56 

0.08 

0.28 

0.08 

0 

Cathode  inlet  gas 

air 

CO, 

Homogeneous 

0.70 

0.30 

rcco 

a(V) 

v*(0)  (V) 

V^(0) (V) 

K 

r  (ft  cm2) 

650 

0.18 

1.03 

1.06 

0.5 

1.08  (PSI) 

1.195  (exp.) 
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f^  =  4-(V*(0)-aa-V*„),  (0<*<l) 
dr  ri,„ 

k(0)=0,  h(  1)  =u, 


where  a  (defined  by  Eq.  ( 8) )  and  rim  are  positive  constants. 

Integration  of  the  differential  equation  with  respect  to  x 
from  zero  to  a  point  in  the  interval  [0,1  ]  gives 


«(X)  = 


v*(Q)-v; 


■M^)) 


(10) 


This  function  must  satisfy  the  end  condition  u(  1 )  —  n,.  Using 
this  condition,  the  following  relation  between  V *cc„  and  u, 
can  be  derived 


v*„=v*(0)- 


Q{1“eKfe)) 


(11) 


The  set  (u(jc),  V^.„)  given  by  Eqs.  (10)  and  (11)  is  the 
solution  of  the  Analytical  Cell  Model.  The  local  current  den¬ 
sity  inside  the  cell  can  be  calculated  from  this  analytical 
solution,  using  Eq.  (6): 


t(x) 


;  dn  V*(0)-V?cn 
‘"6x  r 


(0<x<\) 


(12) 


x 

Fig.  3.  Current  distribution  inside  an  MCFC:  according  to  (cvrve  A)  (Ik 
'analytical  cell  model',  and  ( curve  B)  numerical  computations  based  on  the 
n on-lincarizcd  Nemst  equation  are  shown.  Also  the  constant  (zero-order) 
and  linearly  decreasing  (first-order)  approximation. 

i  .2  **> 

/g)  |(V*(0)-a«-VS„)d« 

0  10  0 


Note  that  Eq.  (11)  is  linear  with  respect  to  the  total  utili¬ 
zation  Thus,  the  use  of  a  linearized  Nernst  equation  results 
in  an  expression  for  the  cell  voltage  that  shows  a  linear 
dependence  on  «].  In  the  following  an  approximation  for 
V£„  will  be  derived  that  is  linear  with  respect  to  the  total 
utilization  «,  as  well  as  to  the  average  current  density  /jnu,. 
We  start  the  derivation  of  this  double  linear  expression  for 
V£.n  with  the  introduction  of  two  more  relations,  i.e.  a  relation 
between  i(0)  and  i ( 1 )  and  a  relation  between  the  cell  voltage 
and  the  electrical  losses  due  to  the  quasi-ohmic  resistance,  r. 

According  to  Eq.  (12)  the  following  identity  holds 

The  cell  voltage  Vfcn  can  be  eliminated  from  this  identity, 
using  Eq.  ( 1 1 ).  Then  the  following  simple  relation  between 
the  current  density  i(0)  at  the  cell  inlet  and  the  current  density 
{( 1 )  at  the  cell  outlet  is  obtained 

aw. 

i(l)-«\0)= - (13) 

r 

In  order  to  obtain  a  relation  between  the  cell  voltage  and 
the  quasi-ohmic  losses  we  multiply  both  sides  of  the  differ¬ 
ential  equation  in  the  analytical  cell  model  by  dw/dx  and 
integrate  with  respect  to  x  from  zero  to  unity 


Rewriting  this  result,  using  the  first  part  of  Eq.  ( 12)  and  the 
following  definition  of  the  average  current  density  <ou, 

(14) 

we  obtain 

i 

f;2(x)d.t  (15) 

2  ,mJ 

Note  that  Eq.  ( 15)  is  exact  within  the  analytical  cell  model. 
This  equation  states  that  the  cell  voltage  is  equal  to  the  aver¬ 
age  value  of  the  linearized  equilibrium  potential  inside  the 
cell,  minus  the  total  loss  in  cell  voltage  due  to  the  quasi-ohmic 
resistance  r. 

It  is  possible  to  calculate  the  integral  in  Eq.  ( 15)  by  using 
the  exponential  relation  for  the  local  current  density  i(x) 
given  by  Eq.  (12).  However,  this  is  not  very  practical  since 
the  expression  for  the  cell  voltage,  obtained  in  this  way,  will 
be  equivalent  to  the  expression  for  V£u  given  in  Eq.  (11). 
Although  the  exact  current  distribution  is  known  to  be  expo¬ 
nential,  it  is  convenient  to  approximate  the  current  density 
i(x)  in  the  last  term  of  Eq.  ( 15)  in  order  to  obtain  approxi¬ 
mate  expressions  for  the  cell  voltage  that  are  more  simple 
than  the  exact  expression  in  Eq.  ( H  ) .  To  this  end,  we  assume 
i(x)  either  to  be  constant  (zero-order  approximation)  or  to 
be  linearly  decreasing  (first-order  approximation).  The 
approximate  current  distributions  must  have  the  same  average 
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current  density  roul  as  the  exact  (i.e.  the  exponential )  current 
distribution  (see  Fig.  3). 

Note  that  the  integral  in  Eq.  (15)  is  equal  to  the  average 
value  of  i2  over  the  whole  cell,  which  is  always  greater  than 
or  equal  to  the  square  of  the  average  current  density,  i.e.  if*,, 
( from  the  Lemma  of  Cauchy-Schwarz  [9] ).  Tills  means  that 
the  last  term  in  Eq.  (15),  that  represents  the  total  loss  in  cell 
voltage  due  to  the  quasi-ohmie  resistance,  r,  is  minimum  for 
a  homogeneous  current  distribution 

i 

—  fr(x)  drkW„ul  (16) 


Z=r— ->0  (19) 

2r/in 

The  significance  of  the  number  Z  will  be  explained  in 
Section  6. 

Comparing  Eqs.  (17)  and  (18)  we  find  that  the  extra 
overpotential  A  i  due  to  the  non-homogeneity  of  the  current 
distribution  inside  the  cell  is  given  by 

A  =~Z2iiout~A2, 

Az  =  0  in  first-order  approximation  ( 20) 


where  the  equality  holds  if  (and  only  if)  the  current  density 
i(x)  is  equal  to  the  average  current  density  rou,  for  all  x.  An 
upper  hound  for  the  cell  voltage  Vj,„  is  thus  obtained  if  we 
assume  a  homogeneous  current  distribution,  i(x)=iiw,  in 
order  to  approximate  the  integral  in  Eq.  ( 15) .  Then  the  zero- 
order  approximation  for  the  cell  voltage  V£n  is  obtained.  This 
approximation  does  not  account  for  the  overpotential,  say  A , , 
due  to  the  non-homogeneity  of  the  current  distribution 

V\*di  =  V *  (0)  - ^au,  - rioul -A„ 


A,  =0  in  zero-order  approximation 


07) 


where  At  >0  can  be  omitted  in  order  to  obtain  the  double 
linear,  zero-order  approximation  for  V*„.  This  double  linear 
relation  for  the  cell  voltage  was  also  derived  by  Machielse 
[  1 J  and  Hcmmes  [  2] ,  using,  however,  different  assumptions 
and  approximations.  The  second  term  on  the  right-hand  side 
of  Eq.  (17)  is  called  the  Nernst  loss  and  is  the  direct  conse¬ 
quence  of  the  decrease  in  V£,(n)  in  the  direction  of  the  cell 
outlet.  The  term  riou,  represents  the  loss  in  cell  voltage  due  to 
the  quasi-ohmic  resistance,  r,  in  the  case  of  a  homogeneous 
current  distribution.  Furthermore,  we  have  shown  that  Aj  is 
the  extra  overpotential  due  to  the  non-homogeneity  of  the 
current  distribution  inside  the  cell. 

Another  possibility  is  the  use  of  a  linear  or  first-order 
approximation  for  the  local  current  density.  According  to  the 
identity  in  Eq.  (13)  this  means  that  »(*)  is  approximated  by 


*(■*)  ~ 


mr, 

r 


first-order  approximation 


Note  that  in  this  first-order  approximation,  the  average 
current  density  is  still  equal  to  ttnil.  Using  the  first-order 
approximation  for  the  local  current  density  to  approximate 
the  integral  in  Eq.  ( 15),  we  obtain  the  following  expression 
forV?cl) 


v*„=v*(0  >--««, -ri„, 


-~Z2riM+A2. 


A2  =  0  in  first-order  approximation  (18) 

where  A2  is  the  inaccuracy  tn  the  first-order  approximation 
and  the  dimensionless  number  Z  is  defined  by 


Although  the  first-order  approximation  is  not  exact,  it  will 
turn  out  that  the  expression  in  Eq.  (20)  with  A2  =  0  is  very 
close  to  the  exact  value  of  A , . 

Furthermore,  it  will  be  illustrated  by  an  example  that  the 
obtained  approximations  for  are  accurate,  notwithstand¬ 
ing  the  fact  that  simple  estimations  for  the  local  current  den¬ 
sity  were  used  in  the  derivations.  Note  that,  once  V£u  has 
been  calculated  using  one  of  the  Eqs.  ( 17),  (18)  or  (II), 
based  on  the  zero-  or  first-order  approximation  for  the  local 
current  density  or  the  exact  solution,  respectively,  a  more 
accurate  expression  for  the  local  current  density  can  be 
obtained  from  Eq.  ( 12). 

We  consider  the  MCFC  as  an  example.  For  an  MCFC 
typical  values  for  V£,(0),  a  and  r  are  V^,(0)  =  1.03  V, 
a  —  1 80  mV  and  r  =  1 .08  Cl  cm2,  respectively.  The  equivalent 
input  current  iin  is  assumed  to  be  200  mA/cm2.  At  a  total 
utilization  of  u,  =  0.75,  the  average  current  density  is  equal 
to  150  mA/cm2,  which  is  a  standard  value  for  an  MCFC 
operation.  Within  the  analytical  cell  model  the  cell  voltage 
V*„  can  be  calculated  exactly  from  Eq.  ( 11 )  to  be  0.7912  V. 
The  zero-order  approximation  in  Eq.  (17)  gives  an  upper 
bound  which  is  9.3  mV  larger,  i.e.  V^„  =  0.8005  -At  V.  This 
means  that  we  have  an  extra  decrease  in  cell  voltage  of  9.3 
mV  due  to  a  non-homogeneous  current  distribution.  The  first- 
order  approximations  in  Eqs.  (18)  and  (20)  turn  out  to  be 
very  accurate.  Eq.  (18)  with  A2  =  0  yields  a  cell  voltage  of 
0.7911  V  compared  with  0.7912  V  for  the  exact  solution. 
Hence,  A2  is,  in  fact,  in  the  order  of  0. 1  mV,  which  is  negli¬ 
gibly  small. 

Fig.  3  shows  the  current  distribution,  obtained  from  Eq. 
(12)  as  well  as  from  numerical  computations  based  on  the 
non-linearized  Nernst  equation.  As  a  consequence  of  the  ini¬ 
tial  dip  in  the  Nernst  potential  (see  Fig.  2  for  small  values  of 
the  local  fuel  utilization)  the  analytically  calculated  current 
distribution  shows  a  much  smaller  current  density  near  the 
cell  inlet,  compared  with  the  numerically  calculated  current 
distribution.  This  implies  that  the  analytically  calculated  cur¬ 
rent  density  will  be  somewhat  larger  in  other  parts  of  the  cell, 
since  the  average  current  density  was  set  equal  for  both 
current  distributions. 

For  the  MCFC  example  (w,=0.75),  the  extra  losses  in 
cell  voltage  due  to  the  non-homogeneity  of  the  current  dis¬ 
tribution  are  small  (in  the  order  of  10  mV).  More  general, 


F.  Standaert  el  at  /Journal  of  Power  Sources  63  (1996)  22J-234 


227 


but  still  within  the  assumptions  of  the  ‘basic  model’,  we  can 
conclude  that  at  normal  (high  efficient)  operation,  the  extra 
electrical  losses  due  to  a  non-homogeneous  current  distribu¬ 
tion  inside  an  isothermal  fuel  cell  are  small  and  will  be  close 
to  the  estimation  (Z2/3)rioul  given  by  Eq.  (20). 

4. 2.  Extended  analytical  cell  model 

In  this  section  we  will  consider  the  influence  of  the  initial 
dip  in  the  Nernst  equation  on  the  cell  voltage.  To  this  end, 
we  will  use  a  partly  linearized  equilibrium  potential,  instead 
of  V*(m).  In  the  so-called  ‘extended  analytical  cell  model’ 
we  consider  the  local  equilibrium  potential  to  be  given  by 
V^u)  for  «  <  u*  and  V^( u)  for  u  >  u*,  i.e.  the  function  that 
follows  the  exact  Nernst  equation  up  to  the  point  of  inflection 
u*  (see  Eq.  ( 8a) )  and  then  connects  with  the  linear  approx¬ 
imation  V*  («)  as  used  in  the  analytical  cell  model. 

It  is  shown  in  Appendix  B  that  the  cell  voltage  cor¬ 
responding  with  this  partly  linearized  Nernst  equation  can  be 
approximated  by 

V$lh  «  V  *  ( 0)  -  -au ,  -  riuul  +  C— , 

2  u  \ 

Mi^O,  zero-order  approximation  (21) 

VS.  =  V*(0)  +|z2j»„lu  +  c-j-, 

U|  =£0,  first-order  approximation  (22) 

with  the  constant  C  equal  to  the  surface  between  the  linear 
approximation  V£,(k)  and  the  initial  dip  in  V^fu) .  Again  the 
zero-  and  first-order  approximations  refer,  respectively,  to  a 
constant  and  linearly  decreasing  current  density  i(x).  Com¬ 
paring  the  approximations  for  and  VSn,  we  see  that  the 
influence  of  the  initial  dip  in  the  Nernst  equation  on  the  cell 
voltage  can  be  taken  into  account  by  an  extra  term,  i.e.  Cl  ut . 

In  the  following  example  it  is  made  clear  that  V**},  is  a  very 
accurate  approximation  for  the  cell  voltage  the  cell 

voltage  corresponding  with  the  non-linearized  Nernst  equa¬ 
tion.  At  least,  as  shown  in  Appendix  B,  for  values  of  «,  in 
the  interval  for  which  the  linear  fit  of  the  equilibrium  potential 
is  accurate  (see  Section  3) 

Vea-VSk  0.2  <«,  <0.8 

The  accuracy  of  both  approximations  in  Eqs.  (21)  and 
(22)  will  be  illustrated  on  the  basis  of  a  similar  MCFC  exam¬ 
ple  as  was  introduced  at  the  end  of  Section  4.1.  However,  we 
now  consider  different  values  for  iin  and  ,  keeping  the  aver¬ 
age  current  density  f0«u  =  *’mMi  constant  at  150  mA/cm2.  For 
the  fuel  composition  as  presented  in  Table  2,  the  value  of  C 
was  numerically  determined  to  be  C  =  3  mV.  Fig.  4  shows 
the  cell  voltage  as  a  function  of  the  total  fuel  utilization,  u,. 
for  a  constant  average  current  density  !«,,=  150  mA/cm2, 
according  to  the  zero-order  approximation  in  Eq.  (21)  as 
well  as  to  the  first-order  approximation  in  Eq.  (22).  The 


Fig.  4.  Example  of  the  cell  voltage  of  an  MCFC  as  a  function  of  the  total 
fuel  utilization,  u,.  for  a  constant  average  cunent  density  of  ISO  mA/cm2 
Top,  middle  and  bottom  curves  correspond  with  the  zero-order  approxima¬ 
tion  in  Eq.  (21 ),  the  first-order  approximation  in  Eq.  (22)  and  numerical 
compulations,  respectively. 

difference  between  these  two  approximations  represents  the 
extra  loss  in  cell  voltage  due  to  the  non-homogeneity  of  the 
current  distribution. 

The  results  of  numerical  computations  based  on  the  non- 
linearized  Nernst  equation  are  also  plotted  in  Fig.  4.  The 
approximations  appear  to  be  accurate  up  to  a  total  fuel  utili¬ 
zation  of  80%.  For  values  of  u,  larger  than  0.8  the  lineariza¬ 
tion  of  the  Nernst  equation  fails  and  the  approximations  start 
deviating. 


5.  Comparison  of  the  models  with  measurements 

In  the  preceding  sections  a  ‘basic  model’  suitable  for  a 
numerical  approach  and  an  ‘analytical  cell  model’  have  been 
introduced.  In  this  section,  some  results  obtained  from  these 
models  will  be  compared  with  experimental  results  presented 
by  Machielse  [  1  ] .  This  comparison  is  mainly  directed  on  the 
cell  voltage  of  an  MCFC  as  a  function  of  the  fuel  utilization, 
in  particular  on  the  following  identity  derived  from  the 
Eqs.  (22),  (14)  and  (19) 


where  the  index  indicates  that  the  partial  derivative  holds 
for  constant  values  of  the  average  current  density 
(*oui  =  constant).  It  will  be  made  clear  that  this  derivative  is 
almost  constant  for  nearly  all  relevant  values  of  m,.  It  strongly 
determines  the  relation  between  the  cell  voltage  and  the  total 
fuel  utilization  ut. 

The  conditions  outlined  in  Section  3  are  the  same  as  used 
by  Machielse.  The  values  a  =180  mV,  C=3  mV  and 
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total  fuel  utilization  Uj 

Fig.  5.  Cell  voltage  as  a  function  of  the  total  fuel  utilization,  U|.  as  obtained 
from  experimental  results  (marked  points)  on  an  MCFC  bench-scale  cell, 
imposing  constant  average  current  densities:  (A)  50  mA/cm2;  (B)  100 

mA/cm2,  and  (C)  150  mA/cm2.  { - )  Plots  for  values  of  h,^0.2  only, 

represent  analytical  results  obtained  from  Eqs.  (21)  (top)  and  (22)  (bot¬ 
tom).  ( - )  Numerical  results  are  plotted  for  all  values  of  u(  <0.9. 

r=  1 .195  H  cm2  will  thus  be  used  in  the  calculations.  For  an 
average  current  density  of  150  mA/cm2  the  three  terms  on 
the  right-hand  side  of  Eq.  (23)  are  successively  —90,  —  15 
and  - 12  mV  for  k,  - 1/2,  the  first  term  being  dominant. 
(Note  that  this  term  is  independent  of  n,,  whereas  the  two 
small  terms  in  the  expression  for  the  derivative  are  not.) 
Nevertheless,  the  sum  of  the  terms  appears  to  be  approxi¬ 
mately  constant  for  a  wide  range  of  values  for  as  is  clear 
from  the  almost  linear  curves  in  Fig.  4.  Calculating  the  aver¬ 
age  value  of  the  derivative  in  Eq.  ( 23)  over  the  interval  from 
«i  =  0.2  to  «i  =  0.8,  we  find 

=  -  124  mV,  analytical;  iwl,5=  150  mA/cm2 

This  value  will  be  compared  with  experimental  data  obtained 
by  Machielse  [  1 1 .  His  experiments  were  carried  out  on  an 
isothermal  1000  cm2  MCFC  bench-scale  cell. 

Machielse  determined  the  cell  voltage  as  a  function  of  the 
total  fuel  utilization,  imposing  constant  average  current  den¬ 
sities.  The  results  are  reproduced,  with  permission,  in  Fig.  5. 
This  figure  shows  three  sets  of  experimental  results,  each  for 
a  constant  average  current  density.  As  can  be  seen,  the  exper¬ 
imental  results  within  a  set  can  be  connected  by  an  almost 
straight  line,  from  which  we  could  determine 

«  —  125  to  —  175  mV,  experimental 

where  the  exact  value  of  the  derivative  depends  on  the 
imposed  value  of  iou1.  Unfortunately,  these  values  are  not 
given  in  Ref.  [1  ].  However,  the  vertical  distances  in  Fig.  5 
between  two  sets  of  measurements  are  proportional  to  the 


current  differences,  because  of  the  term  riimt  in  the  analytical 
approximations  in  Eqs.  (21)  and  (22)  (dashed  curves). 
Thus,  the  average  current  density  /oul  varies  considerably  for 
the  three  sets  of  measurements  and  in  accordance  with  the 
analytical  relation  in  Eq.  (23),  the  slope  of  a  line  through  a 
set  of  experimental  points  depends  on  the  imposed  value  of 
As  the  bottom  set  of  experimental  points  in  Fig.  5  yields 
a  slope  of  approximately  -  125  mV,  very  close  to  the  ana¬ 
lytical  value,  it  is  most  likely  that  this  set  corresponds  to  an 
average  current  density  iout  =  150  mA/cm2. 

For  a  more  detailed  comparison  of  the  experimental  results 
obtained  by  Machielse  and  the  theoretical  results  obtained 
from  either  the  basic  model  ( numerical  computations)  or  the 
extended  analytical  cell  model,  we  assumed  that  the  meas¬ 
urements  were  carried  out  for  average  current  densities  of  50, 
100  and  150  mA/cm2.  For  these  current  densities  the  corre¬ 
sponding  analytical  and  numerical  results  were  presented  in 
Fig.  5  together  with  the  expen  mental  results.  As  shown  by 
this  figure,  the  experimental  and  theoretical  results  are  in  very 
good  agreement. 


5,  Derivation  of  input-output  relations 

As  already  remarked  in  the  introduction,  an  analytical 
approach  can  contribute  to  a  better  understanding  of  fuel  cell 
behaviour.  To  this  end,  some  input-output  relations  will  be 
derived  in  this  section,  i.e.  fuel  cell  performance  will  be 
expressed  in  terms  of  independent  variables. 

Since  the  constant  of  proportion  a  in  the  linearized  Nemst 
equation  can  be  determined  from  the  fuel  composition  (see 
Section  3),  only  two  independent  variables  appear  in  the 
mathematical  formulation  of  the  analytical  cell  model,  i.e.  iin 
and  m,  ( see  Section  4.1).  However,  in  practice,  the  two  ‘con¬ 
trol  knobs’  or  independent  variables  that  govern  the  perform¬ 
ance  of  an  isothermal  fuel  cell  are  the  external  load  resistance 
and  the  equivalent  input  current  density  rin.  Recall  that 
tin  can  be  calculated  from  Faraday’s  law.  For  smaller  test  cells 
one  can  connect  an  electronic  load  that  keeps  the  cell  voltage 
or  total  cell  current  at  a  constant  value  (potentiostatic  or 
galvanostatic  control,  respectively). 

When  the  pair  (/?load,  i,„)  is  considered  as  the  set  of  inde¬ 
pendent  variables,  instead  of  (u|t  iin) ,  it  is  convenient  to  write 
the  cell  voltage  as  the  product  of  the  load  resistance  and  the 
cell  current /oul.  Relating  the  load  resistance  to  the  cell  surface 
BL  we  then  obtain 

V^Rl0iJ0M  =  BLR^i0al  (24) 

On  the  other  hand,  the  cell  voltage  V£.,b  in  the  zero-order 
approximation,  is  given  by  Eq.  (17).  Eliminating  the  cell 
voltage  we  find,  upon  rearrangement 

V*(O)/(flZJiM  +  r+;0 

zero-order  approximation  ( 25 ) 
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This  current  input-output  relation  gives  the  average  current 
density  rtlul  as  a  function  of  the  ‘two  control  knobs’,  ftload  and 
tin* 

Eq.  (25)  looks  like  Ohm’s  law  and  contains  three  ‘resis¬ 
tances’,  all  having  the  dimension  II  m2:  the  scaled  external 
load  resistance,  the  internal  ceil  resistance,  r,  and  a  term  a! 
(2tin)  that  we  call  ‘the  utilization  resistance’.  The  utilization 
resistance  is  due  to  the  Ncrnst  loss  and  decreases  as  the  fuel 
flow  increases.  ( It  vanishes  in  the  limiting  case  of  an  infinitely 
large  fuel  flow,  in  accordance  with  the  fact  that  in  that  case 
fuel  utilization  is  effectively  zero;  so  achieving,  for  a  fixed 
/?)oaa,  the  maximum  output  current  (The  fuel  composi¬ 
tion  is  the  same  throughout  the  cell  and  a  homogeneous  cur¬ 
rent  distribution  is  established.)  In  this  limiting  case  the 
zero-order  approximation  is  exact  and  yields 

W=V*(0)/(5ZJ?load  +  r) 

Recalling  Eq.  ( 19)  wc  see  that  the  dimensionless  number 
Z  was  actually  defined  as  the  rafio  of  two  resistances,  i.e.  the 
utilization  resistance  <*/ (2tin)  and  the  quasi-ohmic  resistance, 
r.  Below  it  will  be  made  clear  that  it  is  the  number  Z  that 
largely  determines  input-output  relations  of  fuel  cells. 

To  obtain  the  output  current  density  as  a  function  of  the 
cell  voltage  and  the  number  Z  we  solve  the  zero-order  approx¬ 
imation  in  Eq.  (  17)  for  u,  and  multiply  the  obtained  equation 
by  iin 


<ou,~  riia+  (a/2)  in 


V£(Q)  -  V&, 
(l+z  )r  ' 


zero-order  approximation 


(26) 


Multiplying  Eq.  (26)  by  V£M  we  obtain  the  power  output 
PoUt  =  per  unit  cell  surface 

V*(0)  -  V*n 

Paul  ~  ctl  — •  LL  V  *el),  zero-order  approximation 

( 1  +  Z)r 


(27) 


Note  that  pou,  as  a  function  of  Vgcl,  is  a  parabola  with  its 
maximum  at  V*^  =  V^(0)/2.  In  Fig.  6  the  parabolas  forZ = 0 
and  Z— 0.5  (i.e.  for  constantly)  are  shown.  The  output  power 
at  constant  utilization  is  given  by  the  equation  that  is  obtained 
if  the  following  identity  is  substituted  into  Eq.  (27) 


This  latter  identity  easily  follows  from  Eq.  (26),  using 
i«  ui =  iin“i =  (2Zr) . 

Maximum  power  output  means  minimum  capital  costs  for 
fuel  cell  stacks.  Consequently,  one  may  be  tempted  to  operate 
a  fuel  cell  at  a  cell  voltage  V£n  =  V£,(0)/2  and  a  fuel  utili¬ 
zation  corresponding  with  the  maximum  of  a  parabola  in 
Fig.  6.  However  from  Eqs.  (24)  and  (25)  it  can  be  shown, 
by  elimination  of  ioul,  that  at  this  point  BLR\mii=r+  (a/ 
2i,n)  =  ( 1  +  Z)r.  Hence,  the  power  loss  due  to  internal  heat 


Veell  (Volt) 

Fig.  6.  Power  output  as  a  function  of  the  cell  voltage  for  constant  values  of 
iln,  hence  of  Z,  as  well  as  for  constant  values  of  u,.  The  circle  denotes  the 
neighborhood  of  normal  operation;  for  parameter  values  see  Table  2. 

production  by  the  quasi-ohmic  resistance,  r,  is  1/(1  +Z) 
times  the  delivered  power  in  the  form  of  electricity,  which  is 
undesirably  high.  Typical  values  for  an  MCFC  at  normal 
(high  efficient)  operation  are  Z~  0.5  and  V^a  =  0.8  V.  Then 
it  follows  from  Eq.  (27)  that  the  output  power  density  pol5, 
can  increase  by  50%  forZ  -»  0,  i.e.  for  tin  -*  w,  even  at  constant 
cell  voltage.  A  further  increase  in  power  output  can  be 
achieved  by  decreasing  the  cell  voltage  by  lowering 
Fig.  6  shows  that  for  small  Z  and  for  a  cell  voltage  equal  to 
half  the  value  of  the  OCV,  a  total  increase  of  100%,  respect 
to  typical  operation,  in  power  output  occurs  (i.e.  from  0.11 
to  0.2-0.25  W/cm2).  However,  at  these  conditions 
BLRfota  ~  r,  which  is  clearly  not  an  optimum  situation  for  fuel 
cell  operation. 

More  accurate  input-output  relations  can  be  obtained  in  a 
similar  way,  using  the  first-order  approximation  in  Eq.  (22). 
For  instance,  analogous  to  the  derivation  of  Eq.  (25)  it  can 
be  shown  that 

+5z2)r+^9] 

4-  —  u i  #  0;  first-order  approximation  (28) 

*  cq(tt) 

Below  we  will  successively  consider  examples  of  the  cur¬ 
rent-input  output  Eqs.  (25)  and  (28),  assuming  a  fixed  load 
resistance.  Examples  of  as  a  function  of  /?!oad  for  a  fixed 
equivalent  input  current  ixn  are  given  in  Ref.  [  10] . 

In  Fig.  7  the  delivered  current  of  an  MCFC  as  given  by  Eq. 
(25),  is  plotted  as  a  function  of  the  variable  for  a  constant 
( and  small)  load  resistace  BLR^^ = 4  ft  cm2.  Again  the  same 
parameter  values,  listed  in  Table  2,  have  been  used. 

Note  that  in  Fig.  7,  the  slope  of  a  line  connecting  a  point 
on  the  curve  and  the  origin  yields  the  total  utilization 
iin  as  indicated  in  the  figure.  Recall  that  the  maximum  output 
current  density  is  achieved  for  an  infinitely  large  value 
of  iin.  Imagine  that,  starting  from  the  point  for  which  fin — 0.8 
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equivalent  input  current  density  iin  (A/cm2) 

Fig.  7.  Output  current  density  as  a  function  of  the  equivalent  input  current 

density  for  an  MCFC  at  a  constant  load  resistance.  ( - )  representing  the 

zero-order  approximation  in  Eq.  (25).  and  ( - )  according  to  numerical 

computations.  Lines  of  constant  utilization  and  the  maximum  achievable 
current,  based  on  an  OCV  of  V^O)  =  1 .06  V,  are  also  shown. 


Fig.  8.  Gose-upof  the  output  current  density  as  a  function  of  the  equivalent 
input  current  density  for  an  MCFC,  showing  very  good  agreement  between 

( - )  the  first-order  approximation  in  Eq.  (28),  and  ( - )  numerical 

computations. 


A/cm2,  the  fuel  supply  will  be  decreased.  As  this  means  that 
iin  decreases,  this  will  result  in  an  increase  in  the  total  fuel 
utilization,  «i-  Since  cannot  be  higher  than  unity  the  line 
k,  =  1  represents  a  limit.  However,  for  large  values  of  the 
linearization  of  the  Nemst  equation  fails  and  consequently 
the  analytical  solution  intersects  the  line  «,  =  1  (as  shown  in 
Fig.  7).  Since  this  is  not  possible,  the  analytical  expression 
in  Eq.  (25)  cannot  be  used  for  very  small  values  of  iia. 
Evidently,  when  uy  is  almost  unity  and  we  continue  to 
decrease  iin,  the  output  current  density  must  approximately 
follow  the  line  «,  =  1  down  to  the  origin.  Numerical  calcu¬ 
lations,  based  on  the  non-linearized  Nemst  equation,  show  a 
quite  sharp  transition  from  the  curved  line  according  to  Eq. 
(25)  to  the  line  u,=  1.  For  iq  >0.8  the  analytical  solution 
starts  deviating  from  the  numerical  solution,  but  significant 
deviations  occur  only  for  very  high  values  of  the  total  fuel 
utilization  ut.  The  analytical  expression  in  Eq.  (25)  is  thus 
very  accurate  up  to  high  utilizations. 

The  output  current  density  i„ul  according  to  Eq.  (28 )  nearly 
coincides  with  the  numerical  solution  up  to  very  high  values 
for  the  total  fuel  utilization  u,  and  is  therefore  not  plotted  in 
Fig.  7.  A  close-up  of  the  curve  obtained  from  Eq.  (28)  is 
given  in  Fig.  8.  This  figure  also  gives  results  for  the  value 
BLR^  =  5  ft  cm2  (bottom  curves) ,  which  is  a  typical  value 
for  the  load  resistance  of  an  MCFC  at  normal  operation 
( foul  —  1 50-1 60  mA/cm2) . 

It  is  evident  that  optimum  conditions  for  fuel  cell  operation 
are  near  the  sharp  transition  point  at  which  the  output  current 
density  diverge  from  the  analytical  solution.  Because  other¬ 
wise,  either  fuel  utilization  becomes  too  high  and  the  cell 
voltage  (given  by  BLff]oadioul)  drops  strongly,  or  the  utiliza¬ 
tion  is  too  low  and  then  not  sufficient  hydrogen  is  converted, 
which  in  general  leads  to  a  lower  system  efficiency.  Fig.  7 


also  shows  that  the  decrease  in  iOM  is  small  up  to  very  high 
values  of  the  total  utilization  Up  Hence,  in  principle,  high 
utilizations  can  be  achieved  without  a  large  decrease  in  cell 
voltage  or  output  power  density  ( i.e.  BL/?,oadi2ut) . 


7.  Conclusions 

The  following  conclusions  could  be  drawn  from  the  present 
study. 

1 .  A  simple  stationary  ‘basic  model’  was  developed,  based 
on  the  assumptions  shown  in  Table  1.  The  ‘basic  model’  can 
be  used  to  describe  isothermal  fuel  cell  behaviour,  independ¬ 
ent  of  the  type  of  fuel  cell. 

2.  It  was  shown  that  the  modeling  of  fuel  cells  can  be 
simplified  by  the  introduction  of: 

•  The  ’equivalent  input  current’,  /iB,  being  the  maximum 
amount  of  current  that  is  produced  in  the  case  of  complete 
fuel  utilization.  It  can  be  calculated  from  the  gas  compo¬ 
sition  and  the  gas  flow  rate,  using  Faraday’s  law. 

•  A 1  inear  relation  for  the  local  Nemst  potential  as  a  function 
of  the  fuel  utilization. 

3.  Linearization  of  the  Nemst  potential  in  the  ‘basic  model' 
results  in  a  cell  model  that  can  be  solved  analytically.  For  an 
MCFC  the  results  obtained  from  the  analytical  cell  model  are 
found  to  be  accurate  up  to  high  fuel  utilizations,  when  com¬ 
pared  with  numerical  computations  based  on  the  non-linear¬ 
ized  Nemst  equation.  The  analytical  cell  model  also  agrees 
well  with  measurements  obtained  by  Machielsc  [  1 J  on  an 
isothermal  1000  cm2  MCFC  bench-scale  cell. 

4.  The  influence  of  simplified  current  distributions  in  the 
cell  was  examined:  a  homogeneous  current  density  and  a 
linearly  decreasing  current  density  were  compared  with  a 
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realistic  current  distribution.  For  an  MCFC  the  loss  in  cell 
voltage  corresponding  with  a  homogeneous  current  distri¬ 
bution  is  almost  10  mV  lower  than  the  actual  loss  in  cell 
voltage. 

5.  The  ‘basic  model’  contains  two  independent  variables, 
which  means  essentially  that  a  fuel  cell  has  two  ‘control 
knobs’,  i.e.  the  external  load  resistance  and  the  fuel 
supply  through  the  equivalent  input  current  /,„. 

6.  The  analytical  cell  model  can  be  used  to  express  the  ceil 
current  and  the  cell  voltage  in  terms  of  the  independent  var¬ 
iables.  Eq.  ( 17)  for  the  cell  voltage  clearly  shows  the  over¬ 
potential  due  to  fuel  utilization.  This  so-called  ‘Nernst  loss’, 
is  proportional  to  the  fuel  utilization  and  the  slope  a  of  the 
linearized  Nernst  equation.  Eq.  (25)  for  the  cell  current  is 
analogous  to  Ohm’s  law.  The  influence  of  the  Nernst  loss 
appears  as  an  ‘utilization  resistance’,  being  proportional  to 
the  slope  a  and  inversely  proportional  to  the  equivalent  input 
current,  /in. 

7.  Fuel  cells  can  deliver  large  peak  powers,  of  course  at 
the  expense  of  losing  efficiency.  Compared  with  normal 
( high  efficient)  operation,  the  power  output  of  an  MCFC  can 
be  increased  by  about  50%  by  just  increasing  the  input  flow 
even  at  constant  cell  voltage.  By  also  decreasing  the  load 
resistance  and  consequently  the  cell  voltage,  the  power 
output  can  be  increased  by  about  100%. 

8.  Fuel  cells  can  be  operated  at  high  utilizations  without  a 
large  decrease  in  cell  voltage  or  output  power  density. 


8.  List  of  symbols 

A ,  decrease  in  cell  voltage  due  to  the  non¬ 

homogeneity  of  the  current  distribution,  V 
A 2  inaccuracy  in  the  first-order  approximation,  V 

B  width  of  the  rectangular  cell,  m 

C  surface  between  the  initial  dip  in  V^u)  and  the 

linearized  function  V£,(«),  V 
c  molar  concentration  of  the  gas  mixture,  mol/m3 

F  constant  of  Faraday,  C/mol 

i  local  current  density.  A/m2 

/in  total  current  equivalent  of  the  fuel  supply,  A 

i*in  total  current  equivalent  of  the  fuel  supply  per  unit 

cell  surface.  A/m2 

/ou,  total  current  output  of  a  fuel  cell,  A 
foul  total  current  output  of  a  fuel  cell  per  unit  cell 
surface;  average  current  density.  A/m2 
hm*  maximum  output  cunent  density  for  a  fixed  value 

of  the  load  resistance.  A/m2 

K  equilibrium  constant  for  concentrations  in  the 

shift  reaction 

L  length  of  the  rectangular  cell,  m 

/t  total  number  of  moles  in  a  control  volume,  mol 

n,  number  of  moles  of  a  species  in  a  control 

volume,  mol 

p„u,  output  power  density,  W/m2 


Pi  partial  pressure  of  the  <th  species  in  the  fuel,  N/ 

m2 

r  constant  quasi-ohmic  resistance,  li  mr 

/?louJ  external  load  resistance,  H 

R  universal  gas  constant,  J  /  ( mol  K) 

T  absolute  cell  temperature,  K 

u  fuel  utilization 

u*  point  of  inflection  of  the  function  V^u) 
k,  total  fuel  utilization,  i.e.  value  of  u  at  the  cell 

outlet 

v  utilization  of  CO 

Vggi  u)  Nernst  potential  as  a  function  of  the  utilization  #, 

V 

V£j(w)  linear  fit  for  V„(ii),  V 

cell  voltage,  V 

Vco,  intercept  of  linearized  polarization  curve 

( Kor = 0  for  the  MCFC),  V 

Vfcl,  cell  voltage  if  the  Nernst  equation  is  linear  with 
respect  to  u,  V 

V“u  cell  voltage  corresponding  with  a  partly 

linearized  Nernst  equation,  V 
x  scaled  (i.e.  relative)  distance  to  the  cell  inlet 

X  distance  to  the  cell  inlet,  m 

X,  mole  fraction  of  the  ith  species 

Z  dimensionless  number  (see  Eq.  ( 19) ) 

Greek  letters 

a  characteristic  value  for  the  derivative  of 

—  VcqC")  with  respect  to  «,  V 
tj  local  overpotential,  V 

£  integration  variable 

Abbreviations 

OCV  open-circuit  voltage,  i.e.  V’eq(0) 

MCFC  molten  carbonate  fuel  cell 

ODE  ordinary  differential  equation 

BVP  boundary  value  problem 


Appendix  A,  The  equilibrium  (or  Nernst)  potential  as  a 
function  of  the  utilization  in  the  case  of  an  MCFC 

In  this  appendix  an  expression  for  the  local  equilibrium 
potential  inside  an  MCFC  as  a  function  of  the  local  fuel 
utilization,  u,  is  derived.  The  derivation  is  based  on  a  common 
representation  of  the  Nernst  equation  and  on  the  assumptions 
2,  3,  6  and  7  given  in  Table  1.  Furthermore,  we  neglect 
diffusion  effects  in  the  gas  phase,  i.e.  the  anode  gases  are 
considered  to  be  homogeneous  over  the  cross  section  of  a 
channel  (conform  to  assumption  4).  We  also  neglect  the 
methane  reforming  reaction  and  distinguish  only  four  chem¬ 
ical  species  in  the  fuel  flow,  i.e.  H2,CO,  H20  andC02.  Partial 
properties  are  denoted  with  an  index  is  { 1,...,4J,  using  the 
following  sequence  for  the  species 
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Species  H2  CO  H20  C02 

The  cell  half-reaction  at  the  anode  side  of  an  MCFC  is 
given  by 

H2+C032  ^H2O  +  C02  +  2e-  (Al) 

Assuming  constant  partial  pressures  for  the  oxidant  spe¬ 
cies,  the  corresponding  Nemst  equation  can  be  written  as 


V^iuy^V^Q) 


RT  MO)  T  Pl(u ) 

IF  n^3(0)p4(0)  IF  np3(«)p4(«) 

(A2) 


where  V^O)  is  the  OCV.  The  OCV  can  be  calculated  from 
the  inlet  gas  compositions  which  are  assumed  to  be  known. 
If  the  fuel  behaves  like  an  ideal  gas  we  can  use  the  idea!  gas 
law  to  determine  the  partial  pressures,  p/2  of  the  components 

p  =cRTXh  i  —  1,  ...,4  (A3) 


«,(«)+n2(K)  =  (rt,(0)  +n2(0))(l  ~u)  (A7) 

The  total  number  of  moles  CO  that  is  converted  by  the 
shift  reaction  if  complete  fuel  utilization  is  achieved  is  given 
by  n2(0).  The  fraction  of  this  number  that  is  converted 
between  the  cell  inlet  and  the  position  that  corresponds  with 
the  utilization  u,  is  denoted  by  v(u).  Hence,  by  the  definition 
of  o  we  also  have 

n2(M)=«2(°)(l“y(«))  (A8) 

Below,  the  mole  fractions  X,,  as  given  by  Eq.  (A5),  are 
expressed  in  terms  of  u  and  u(w).  The  equilibrium  potential 
in  terms  of  u  and  u(k)  follows  from  the  representation  of  the 
Nernst  equation  in  Eq.  ( A4) .  Next,  the  equilibrium  condition 
for  the  shift  reaction  is  used  to  determine  v  as  a  function  of 
w,  after  which  the  equilibrium  potential  can  be  expressed  in 
terms  of  u  only. 

In  order  to  determine  the  mole  fractions  as  functions  of  u 
andu(it)  we  first  subtract  Eq.  (A8)  fromEq.  (A7),  yielding 


with  X,  the  mole  fraction  of  the  rth  species,  c  the  molar  con¬ 
centration  of  the  mixture  and  R  the  universal  gas  constant. 
Note  that  c  is  constant  throughout  the  cell,  since  we  assumed 
temperature  and  pressure  to  be  uniform.  Using  Eq.  (A3)  to 
eliminate  the  partial  pressures  from  the  Nernst  Eq.  (A2)  we 
obtain 


Vcti(u)=Vai(Q)~\n- 


X,(0) 


RT, 

-+~—  In- 


X,(«) 


2  F  X3(0)X4(0)  2  F  X3(«)X4(«) 

(A4) 


Hence,  to  express  the  local  equilibrium  potential  in  terms  of 
u  we  have  to  determine  the  mole  fractions  X,  as  functions  of 
u. 

To  this  end,  the  number  «*(()),  ( f  =  1 .... ,4 )  is  defined  as 
the  number  of  moles  of  the  rth  species  that  enters  the  cell  at 
the  anode  side  per  unit  of  time.  The  total  number  of  moles 
n(0)  that  enters  at  the  anode  side,  per  unit  of  time,  is  given 
by  the  sum  of  «,(0)  over  all  i.  The  numbers  «,(«)  and  n(«) 
are  defined  in  an  analogous  way  and  are  relatea  to  one  and 
the  same  volume  element,  following  the  motion  of  the  flow. 
The  mole  fractions  of  the  species  as  a  function  of  the  utili¬ 
zation  can  be  written  as 


1-1 . 4  ( A5) 

n(n) 

The  numbers  «<(«)  in  Eq.  (A5)  depend  not  only  on  the 
conversion  of  hydrogen  by  the  electrochemical  reaction,  but 
also  on  the  production  of  hydrogen  due  to  the  water-gas  shift 
reaction 

CO  +  H2Os*C02  +  H2  (A6) 

The  total  number  of  moles  H2  that  can  be  converted  by  the 
electrochemical  reaction  in  the  case  of  complete  fuel  utili¬ 
zation  is  equal  to  the  number  of  moles  H2  and  CO  that  enters 
the  cell.  Since  the  utilization  n(x)  is  defined  as  the  fraction 
of  this  supply  that  is  converted  in  the  interval  [0*xl  we  have 


/!,(«)  =n,(0)  -  (n,(0)  +n2(0))K  +  n2(0)D(u)  (A9) 

The  second  term  in  the  right-hand  side  of  Eq.  (A9)  gives 
the  decrease  in  the  number  h,(m)  due  to  the  electrochemical 
reaction.  Hence,  from  the  reaction  Eq.  (Al)  it  follows  that 
in  the  control  volume  (n,(0)  +n2(0)  )u  moles  C02  and  as 
many  moles  H20  are  produced  by  the  electrochemical  reac¬ 
tion.  The  last  term  in  Eq.  (A9)  is  equal  to  the  number  of 
moles  H2  produced  by  the  shift  reaction.  From  tne  reaction 
in  Eq,  ( A6)  it  follows  that  in  the  control  volume  n2(0)u(u) 
moles  C02  are  produced  and  as  many  moles  H20  are  con¬ 
verted  by  the  shift  reaction.  Hence,  the  numbers  «3(  n)  and 
«4(m)  are  given  by,  respectively 

n3(u)  =n3(0)  +  (n,(Q)  +n2(0)  )w-/i2(0)u(w) 

iz4(«)  =.n4(0)  +  (n,(  0)  +n2(0))M+/i2(0)u(«)  (A10) 

The  U  tal  number  of  moles  in  the  control  volume  can  now 
be  calculated  by  adding  the  numbers  n, 

4 

n(u)  =  £*,(«)  =n(0)(l  +  [X((0)  +X2(0))«)  (All ) 


Using  the  identities  in  Eqs.  (A5),  (A7)-(A11)  and  the 
following  notation 


X,(0)  X,(0)+X2(0) 

X2(0)-  X2(0) 


1  +at. 


<Mn)  =  l  +  (X,(0)+X2(0 ))«  ( A12) 

The  mole  fractions  in  Eq.  ( A5)  can  be  written  as 


X,<«)  =  X,(0)^1  — u)j!  <M« 
X,(ii)  —  X2(0)(1  —  u)/<£(«) 

X,(«>  =x,(0)|l  +!<*-.,)/, 


<M  ») 
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x4(«)  -X,(0)(l  h1'))/ +00  (A13) 

We  assume  (hat  the  water-gas  shift  reaction  in  the  MCFC 
is  in  equilibrium  throughout  the  cell.  Using  the  identities  in 
Eq.  (A3)  orEq.  (A  12),  the  equilibrium  constant  K  for  to  the 
reverse  shift  reaction  can  be  rewritten  as,  respectively 

Pi>P3_Xa(0)X>(0)  X2(u)X,(u) 

P  \P4  X,(0)X4(0)  X,(a)X4(«) 
or 


K=- 


Substitution  of  the  expressions  for  the  mole  fractions 
(Eq.  (A  13))  into  the  equilibrium  condition  yields 

=  (l-^(*K-ll))(l+^-(i>K  +  L0)  <AI4) 

Solving  Eq.  (A14)  for  v(u)  and  using  the  equilibrium  con¬ 
dition  again,  we  find 

b 


u(u)  =d  + 


2(1  ~K) 


(.  i  Y  V  .  (2 Ka,  +  a,)6 

rr  johoW  +—~u- 


(0<K<1)  (A15) 

where  the  trivial  condition  u(0)=0  is  satisfied  and  d  is  a 
constant  defined  by 


d= 


1  +a$  +  K(a\  +a4) 
2(1 -if) 


Using  Eqs.  (A13)  the  Nemst  Eq.  (A4),  for  zero  oxidant 
utilization,  can  finally  be  rewritten  in  terms  of  u 


V«,(«)  =  Vcq(0)  + 


X7T 
2F  [ 


j^ln^l  —^-(bu  —  u(n))j+ln  d>(n) 
-ln|l  +  j-(fcu-t>(n))j-ln^l  +^-(&u  +  i/(u))jj 


(A16) 


utilizing  for  $(«)  and  i>(n)  the  expressions  given  by  Eqs. 
(A12)  and  (A15),  respectively. 

In  an  analogous  way,  the  obtained  expression  for  Vcq(«) 
can  be  extended  to  an  expression  that  also  accounts  for  the 
utilization  of  oxidant.  However,  this  latter  step  is  also 
described  in  Ref.  [4]  and  has  been  omitted  in  this  paper. 


A.  I.  Nota  bene 


Fig.  9  shows  an  example  of  the  reaction  coordinate  w(w) 
for  the  conversion  of  CO  by  the  shift  reaction.  For  small 


local  fuel  utilization  u 

Fig.  9.  Reaction  coordinate  for  the  conversion  of  CO  by  the  shift  reaction, 
in  the  case  of  standard  gas  compositions  (see  Table  2). 


values  of  the  fuel  utilization,  CO  is  produced  from  the  reverse 
shift  reaction,  while  when  the  hydrogen  partial  pressure 
becomes  low,  the  shift  reaction  evolves  towards  the  decom¬ 
position  of  the  produced  CO.  The  value  of  n>  0  for  which  v 
is  equal  to  zero  is  obtained  by  setting  u =0  in  Eq.  ( A14)  and 
solving  the  result  for  u 


(1  +  1/X)X2(0)+X4(0) 
X,(0)+X2(0) 


u(u)=0 


This  value  is  exactly  0.5  for  the  standard  gas  compositions 
(see  Table  2)  and  AT =0.5. 


Appendix  B.  Cell  voltage  as  a  function  of  the  totai 
utilization 

In  Eq.  ( 15)  the  cell  voltage  V£„  in  the  solution  of  the 
analytical  cell  model  was  written  as  the  average  value  of  the 
linearized  equilibrium  potential  minus  the  polarization  losses 
due  to  the  quasi-ohmic  resistance,  r.  In  the  same  way  Eq. 
(15)  was  derived,  but  starting  with  the  ODE  in  the  BVP,  it 
can  be  shown  that  the  cell  voltage  Vte11  in  Section  2  can  be 
written  as 

«■  i 

Kc, = -  f IV «)  du  - — f r ( *)  dt  (AI7) 

J  foul* 

0  0 

In  Section  4.1,  it  was  shown  (see  Fig.  3)  that  the  initial 
dip  in  the  Nemst  equation  has  a  very  small  influence  on  the 
current  distribution  inside  the  cell.  Hence,  at  normal  opera¬ 
tion,  the  second  integral  in  Eq.  (A17)  can  still  be  approxi¬ 
mated  using  a  constant  (zero-order)  or  linear  (first-order) 
approximation  for  the  current  distribution  inside  the  cell  Just 
as  in  the  reduction  of  Eq.  ( 15).  Using  the  first-order  approx¬ 
imation,  Eq.  (A17)  can  be  reduced  to 
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first-order  approximation  (  A 18 ) 

Assuming  that  V^( u) ,  see  Section  3,  is  an  accurate  approx¬ 
imation  of  ^(tt)  for  values  of  u  in  the  interval  (u*,w,),  we 
set  Voq(u)  inEq.  (A18)  equal  to  V*(w)  for  values  of  u>u*. 
In  this  way  we  obtain  the  following  expression  for  the  cell 
voltage 

v«„  =  v*  (0)  -(l  +iz2)r/„11,  +  C~,  u,  >  u* 

(AI9) 

where  C  is  the  surface  between  the  initial  dip  in  the  equilib¬ 
rium  potential  and  the  linear  fit  V£,(a) 

c=  f[V,„(«)-V*(ii))dK 

(> 

The  approximation  in  Eq.  (A  19)  is  also  valid  for  values 
of  u i  <  «*,  as  long  as  V£,(  u)  is  an  accurate  approximation  of 
VcqfK)  for  values  of  u  in  the  interval  (u,,u*).  Hence,  Eq. 
(A19)  is  valid  for  values  of  u,  in  the  interval  for  which  the 
linear  approximation  V£,(u)  is  accurate. 
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